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Abstract 



In a class of F-theory SU{5) GUTs the low energy chiral mass spectrum is obtained 
from rank one fermion mass textures with a hierarchical structure organised by U{1) 
symmetries embedded in the exceptional Eg group. In these theories chiral fields reside 
on matter 'curves' and the tree level masses are computed from integrals of overlapping 
wavefuctions of the particles at the triple intersection points. This calculation requires 
knowledge of the exact form of the wavefuctions. In this work we propose a way to 
obtain a reliable estimate of the various quantities which determine the strength of the 
Yukawa couplings. We use previous analysis of KK threshold effects to determine the 
(ratios of) heavy mass scales of the theory which are involved in the normalization of 
the wavefunctions. We consider similar effects from the chiral spectrum of these models 
and discuss possible constraints on the emerging matter content. In this approach, we 
find that the Yukawa couplings can be determined solely from the U{1) charges of the 
states in the 'intersection' and the torsion which is a topological invariant quantity. 
We apply the results to a viable SU{5) model with minimal spectrum which satisfies 
all the constraints imposed by our analysis. We use renormalization group analysis 
to estimate the top and bottom masses and find that they are in agreement with the 
experimental values. 



1 Introduction 



Recent progress in F-theory model building [1 h8J I has shown that old successful GUTs, 
including the Georgi-Glashow minimal SU(5), the SO(IO) model etc, are naturally 
realised on the world-volume of non-perturbative seven branes wrapping appropriate 
compact surfaces. The rather interesting fact in F-theory constructions is that they are 
defined on a compact elliptically fibered Calabi-Yau complex four dimensional manifold 
thus the exceptional groups Eq,Ei,E^, can be naturally incorporated into the theory 
too [IHSllB]. Although exceptional gauge symmetries suffer from several drawbacks 
when realized in the context of four-dimensional grand unified theories, in the case of 
F-theory models they are more promising as new possibilities arise for the symmetry 
breaking mechanisms and the derivation of the desired massless spectrum. 

Present studies on F-theory model building have been concentrated on three gener- 
ation -mainly SU{5)- GUT models which fall into the following two distinct categories: 
those where all three families with the same Standard Model representation content 
are assigned to a single matter curve |llH14j . and variants [15H21] where some or all 
of the quark and lepton families are assigned to different curves. Several of these con- 
structions built up to these days have attempted to give solutions to fundamental GUT 
problems as is the case of doublet-triplet splitting, the rapid proton decay, the Higgs 
mixing term, the neutrino sector and other related issues [3l[6l [T0l[16ll21f[23] . To analyse 
the phenomenological properties one should extract the relevant information from the 
superpotential which can be readily constructed once a particular assignment of the 
fermion families and Higgs on the matter curves has been chosen. Of course, domi- 
nant role on the estimation of such effects is played by the Yukawa couplings, thus the 
theory's predictive power depends on the calculability of the latter. 

In F-theory GUTs the trilinear Yukawa couplings are realised at the intersections 
of three matter curves Sj, i = 1,2,3 where the zero-modes of two fermion fields and a 
Higgs boson reside. Along these curves the Gs symmetry is enhanced Gs^ 3 Gs x U{l)i 
while the corresponding zero modes are charged under the U{l)i. To determine the most 
general structure of the zero-mode wavefunctions one has to solve their corresponding 
differential equations of motion emerging from the twisted eight-dimensional Yang- 
Mills action [2], (see also [14ll24f[26] ). In general, the solutions are found to exhibit the 
expected gaussian |24p25j profile which falls off exponentially away from the curve while 
their exact form is specified by a mass scale characterizing the size of the compact space 
and the particular C/(l)i-charge of the relevant zero-mode. The Yukawa couplings of 
the {33}-entries of the up, down and charged lepton mass matrices are then computed 
in terms of the integrals of overlapping wavefunctions of the aforementioned form at 
the intersection point of three matter curves 



Js 

where the integration is over the compact internal four-dimensional manifold S sup- 
porting the gauge group of the theory, ip, (p denote the fermion and Higgs wavefunctions 

^ For other recent related work and a review see also [91110) 
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respectively whilst M^, is associated to the supergravity limit of the compactified the- 
ory, hence it will be naturally linked to the F-theory compactification scale. Further- 
more, assuming that higher order non-renormalizable Yukawa couplings are generated 
through mediation of heavy string modes and/or Kaluza-Klein states, the calculation 
of all Yukawa entries of this type in the fermion mass matrices can be reduced to a 
similar computation [18]. 

To quantify low energy implications of a given model, a reliable estimation of the 
integral ([T|) is required. This is however a non-trivial task since, although the local 
profile of the wavefunction is precisely known, the final result is expressed in terms of 
several parameters which mainly appear in the normalization of the wavefunction. In 
particular, the background scalar field vev, the Mqut scale, the scale M* as well as 
the geometry of the surface S supporting the GUT symmetry are not precisely known. 
Estimates of such global quantities have been given and rather sensible results have 
been obtained in the case of the mass hierarchy [12lll5t [T8]. However, in this context 
one fails to predict a heavy top mass compatible with the experimental data in an 
unambiguous manner. 

In this work we reconsider the problem of trilinear Yukawa couplings in F-theory 
GUTs, focussing on the normalization of the wavefunction and the role of the various 
scales of the theory. We rely on the results of [27] and [5] to determine the effects 
of KK-massive modes on the renormalization group running the GUT scale and the 
gauge unification. Using these results we estimate the modifications of mass scales 
ratio Mgut/Mc and express this solely in terms of the Ray-Singer torsion [28] which 
is a topologically invariant quantity. Identifying the various scales appearing in the 
computation of the overlapping integrals with Mqut and Mc, we argue that this gives 
a reliable computation for the top Yukawa coupling. The final result for the Yukawa 
coupling depends only on the charges of the fields involved and an exponential factor 
involving the torsion. 

In the next section we review in brief the basic F-theory set up and the relevance 
of the Ray-Singer torsion to the threshold corrections. We show how KK-massive 
modes modify the GUT scale and in section 3 we perform a similar analysis for the 
case of chiral matter subject to flux constraints. In section 4 we apply the analysis 
in a realistic F-theory SU{5) GUT [18j which was shown to reproduce successfully 
the fermion hierarchy and confront other main phenomenological issues. In sections 
5 and 6 we perform the renormalization group analysis and obtain the formulae for 
the top and bottom Yukawa couplings by computing the corresponding integrals. As 
already said the result is expressed in terms of the analytic torsion. As an application, 
we perform an explicit numerical calculation of the top and bottom Yukawa couplings 
choosing a specific line bundle on a Hirzebruch surface. In section 7 we summarize 
our work and present our conclusion. Some technical details regarding the solutions of 
the renormalization group equations for the computation of the third generation quark 
masses are shown in the appendix. 
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2 KK-modes and the GUT scale 



Our F-theory set up is assumed to be the world- volume of the seven-brane of a ADE- 
type singularity which wraps the space R'^'^ x S where is a Kahler manifold of two 
complex dimensions. At low energies, F-theory is described by an eight-dimensional 
Yang-Mills theory on R^'^ x S which must be topologically twisted to preserve = 1 
supersymmetry [2l|3]. The super symmetric multiplet consists of the gauge field, two 
scalars <^8^9 combined into the complex fields = <I>8ibi<I>9 and the fermions r/, -0, x ^ 
in the adjoint of the gauge group. In the twisted theory ry, x appear as holomorphic 
(/c, 0) forms (and their complex conjugate antiholomorphic (0,?)) with k,l = 0,1,2, 
while the scalar (/> is a two form. At the d = A, N = 1 level, these fields are organised 
as one gauge and two chiral multiplets 

(012, XI2) 

For the computation of threshold corrections we assume that the gauge theory 
model is described by some GUT gauge group which for definiteness we take it to be 
Gs = SU{5) (at a scale Mqut ^ Mpianck)- In this limit it is also natural to assume 
that KK-modes are much lighter than other string excitations. 

Threshold corrections are expressed in terms of the masses of the KK-spectrum 
which emerge as the non-zero eigenvalues of the Laplacian of the eight dimensional 
theory. Therefore, threshold corrections are constants and should be independent of the 
particular metric on the surface 5, while they are expected to be expressed in terms of 
a topologically invariant quantity. It was first shown, already sometime ago [27], that 
in the context of M-theory compactified on a manifold with G2 holonomy the most 
appropriate topological invariant quantity to express these threshold corrections is the 
Ray-Singer torsion [28j. A similar analysis has been performed for the case of F-theory 
in j5] which we now review in brief. 

The SU (5) subgroup which commutes with the Standard Model {Gsm) is the U (l)y . 
We denote with qi the U{1)y charge of the Standard Model representations Ri arising 
from the decomposition of the adjoint of SU{5) under Gsm 

24 ^ (8, 1)0 + (1, 3)0 + (1, 1)0 + (3, 2)_5/6 + (3, 2)5/6 (2) 

We assume that the GUT group breaks to Gsm by turning on a flux along U (l)y , hence 
each Qi determines now a line bundle which we denote it by . We observe that in the 
SU{5) adjoint decomposition ([2]), in addition to the standard SU{3), SU{2) adjoints 
there are in principle exotic zero modes in the (3, 2) and (3, 2) representations along 
the £='=5/6 ijj^g bundles whose numbers are given by the Euler character —xiS,C^^^^). 
Therefore, elimination of these massless modes requires x{S, C^^^^) = which implies^ 
the following relation for the cohomology class: ci(£^/^) • ci(£^/^) = — 2. 

Consider next the Dolbault operator d of the corresponding holomorphic bundle V 
with representation R{V), acting on /c- forms as follows 

8 : n°^'' ® R{V) R{V), for = 0, 1 (3) 

^The number of massless modes of a 7i representation is —x{S,Ti) = 1 + |ci(7i) • (ci(7i) + ci(S)) 
while here the following relation is satisfied: ci{£^^^) ■ ci{S) = [S]. 
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and the Laplacian 

A,^n(v) = {d + d^f = dd^ + d^d (4) 
If we denote collectively with ip^ its A;-form eigenfunction then 

AkMV)rk = A^Vr (5) 

where represents the corresponding eigenvalue and in four dimensions corresponds 
to a squared mass. 

We express now the threshold corrections of the gauge multiplet in terms of the 
Ray-Singer torsion. The running of the gauge couplings are governed by the equation 

-^ = —^ + Klog-^+Sjf\ « = 3,2,y (6) 

i/a \H'J its H- 

where A is the gauge theory cutoff scale, ka = (1, 1,5/3) are the normalization coeffi- 
cients for the usual embedding of the Standard Model to SU{5), Qs is the gauge coupling 
as deduced from the higher theory and ha are the one-loop /3-function coefficients which 
are given by H 

6, = 2StrA/=oQl(3^-x') (7) 

where x is helicity operator, Str denotes the supertrace (bosons contribute with 
weight -|-1 and fermions with —1) and Qa stands for the three generators of the Standard 
Model gauge group ^[/(S) X SU{2) X C/(l)y. Finally, 5^^^ correspond to the one-loop 
threshold corrections and are given by a similar term over the massive states (in contrast 
with the previous term that corresponds to massless states) 

= 2 StvM^oQl - X') log(AVM2) (8) 

Ri 

= 2 ^ Trn^Ql StrM^o - x') \og{K^/M^) (9) 

In the second line StrM^oQa been factored out since it depends only on the rep- 
resentation Ri (we are suppressing the notation for the bundle V in RiiV)). For a 
certain helicity state there is a logarithmic dependence on its mass squared which 
corresponds to the eigenvalue of the Laplacian /j,. Further, since the trace of log is 
equal to the logarithm of the determinant we can replace the trace of log A^/M^ above 
with 

-logdet'(Afc,H,/A2) (10) 



where the prime on det denotes that we exclude the zero modes. 

Now, each eigenvector of the zero-form Laplacian Aq r contributes a vector mul- 

2 ' 2 ' 



tiplet with helicities 1,-1, while the one-form Laplacian Ai^/j. gives a chiral 



^We adopt here the notation of [27] , 
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multiplet with helicities 0,0,^,—^. Similarly, ^2,Ri is associated to anti-chiral multi- 
plets. Evaluating the supertrace we get 

/ N r — I for the vectot multiplet 

Str(--x')= (11) 
^ ^ I +^ ^'^^ chiral multiplet 

Therefore, the total sum of the contribution of KK-modes to the thresholds from the 
gauge fields is then written 

5(^)=2^Tr^,(Q2)/c, (12) 



with 



and 



/Cj = -logdet — -V^ logdet — logdet — (13) 



(n2 \ _/ {3,2,0}, for (8,l)o + (l,3)o + (l,l)o 
Tri^o ^Q{3,2,y}j -t{l,3,f}, for (3, 2)_5/6 , (3, 2)^/6 ^'^^ 

with normalization Tr((5^) = Therefore, the three 5^^'* read 



^ = y (^5/6 + ^-5/6 J 
5('^) = 4/Co + 3(K5/6 + /C_5/6) (15) 

5^^) = 6/Co + 2(/C5/6+/C„5/6) 

with an obvious notation in the subscript of /C. Further, the spectrum of Ai is equiv- 
alent to the sum of the spectra of Aq and A2. Then we can write (jlSp as follows 



/Ci = 21ogdet'^J#i-logdet'^i#i (16) 

A^ 

Introducing the Ray-Singer torsion 

n = \j2{-lf^'\ogdei'^ (17) 

fc=0 

and noting that the Laplacian commutes with the Hodge * operator which maps A;-forms 
to (2 — fe)-forms, so that A^^/j and lS.2^k,R* have the same spectrum, we get ^ 

Tr = -Tr* = \ (2 logdet'^ - logdet'^) (18) 

The expression inside the parenthesis is exactly what we have derived previously for 
the threshold corrections. According to the Ray and Singer theorem |28J the torsion is 
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independent of the metric of the manifold and therefore independent of the cutoff scale 
A. (The latter can be eliminated by scaling the metric.) 

In the above expression for the torsion it has been assumed that there are no zero 
modes as is the case of a non-trivial representation i?,. However, in the case of the 
trivial bundle there are zero modes and the torsion is not equal to /C^q. These are 
related by 

2To = lCo-\og{V h^) (19) 

where in the subscript we used the standard notation O for the trivial representation. 

We assume 5 to be a del Pezzo surface and we expect that /C5/6 = ^-5/6 so that 
the equations in ()15p can be cast in the form 



^.^^\s^^\si'^) = (y/Cs/e, 6/C5/6 + 4/Co, 4/C5/6 + 6/Co) (20) 



Taking into account that the contribution of the gauge multiplet to the /3-function is 
b^a"* = (0, —6, —9) we can rewrite the above as 

Si'^ = (/C5/6 - /Co) + 10A;,/C5/6 = (Ts/e - To) + 20 k,%„ (21) 

The second term is proportional to ka and therefore it can be absorbed to the IGvr^/ca/g^ 
term, redefining the coupling gs- And, finally, the initial equation ^ for the running of 
the gauge coupling, where only the gauge multiplet has been taken into account, gives 



levr^fc, ^ exp[4/3(r5/6-ro)] 

= — + log 2T7l72 

where all the A dependence, after the appropriate cancellation between the massless 
and massive modes, has been absorbed in the redefinition of the bare gauge coupling 
along with the second term of (j2ip . and we can define the Mqjjx as 



^2 _ exp [4/3 (7^/6 -7^))] 



M6uT = (23) 

If we associate the world volume factor with the characteristic compactification 

scale Mc, we can write this equation as follows 

Mgut = e^/^('^5/6-7?)) . (24) 

Thereby, we conclude that the ratio Mgut/Mq of the two characteristic scales of the 
theory depends only on the torsion and is topologically invariant. 
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3 The inclusion of chiral matter 



We have seen in the previous section that the contribution of threshold corrections 
from KK-modes of the vector multiplet leads to a cutoff independent RG running of 
the gauge couplings. We wish now to include the analogous contributions from the 
chiral and Higgs sector of the theory. In F-theory constructions we need to take into 
account the zero mode as well as the KK-massive mode contributions. 

We first start by incorporating the zero mode effects. In the case of the SU (5) theory 
we discuss in this paper the possible non-trivial representations transform as 10, 10 and 
5, 5. Matter fields in general, will contribute to the RG running in (j22p with a term 

. ;2 

of the form 6^ log where 5^ are the /3-function coefficients for the corresponding 
three SM gauge group factors and A' a cutoff scale not necessarily equal to the gauge 
cutoff A. If matter fields in the theory arise in complete SU{5) multiplets, then the 
/3-functions contribute in proportion to the coefficients ka, ie. oc ka- Then writing 

log ^ as log h log ^'^■F^ we can absorb the first logarithm in a redefinition of 

the gauge coupling 

^ ka = ka = —ka + bl log — 5 (25) 

9gut o,gut as ^^GUT 



The massless chiral and Higgs spectrum at low energies, however, does not form com- 
plete SU{h) multiplets. As it is well known, Higgs doublets are found in 5 -|- 5 rep- 
resentations of SU{5) together with the dangerous color triplets. To protect proton 
from rapid decay, triplets must receive a mass at a relatively high scale Mx < Mgut- 
Taking this into account we write = ba + b^ where ba denotes the MSSM /^-functions 
while b^ the color triplet part. Adding all contributions, ()22p takes the form 

It should be mentioned that if the contributions from matter fields fail to obey the 
condition b^ oc ka, then ([25]) splits the common gauge coupling to three different 
values at Mgut- 

We discuss now the threshold contributions arising from the Sg and Sio matter 
curves. From the decompositions 

10 ^ (3, 2) 1 + (3, 1)_2 + (1, l)i, 5 ^ (3, 1) 1 + (1, 2)_i 

6 3 3 2 

we readily find that KK-modes residing on these curves contribute to thresholds as 
follows ^ 



Sa = ^-1/2, ^-1/2 + 2/3/C1/3) 

Sa^ = (2/C1/6 + ^-2/3, 3/C1/6, I/3/C1/6 + 2/Ci + 8/3/C_2/3) 



*In the next sections we will see however that gauge coupling splitting at Mqut is not avoided when 
a non-trivial U{1)y flux is turned on to break SU{S). 
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where, as previously, /Cj stand for the Ray-Singer torsion. As in the case of gauge 
fields one can absorb the parts proportional to ka by redefinition of the initial gauge 
coupling. Further, it is anticipated that on appropriate bundle structures such correc- 
tions will diminish or be negligible compared to other (from light degrees of freedom 
or supersymmetry) threshold effects [22]- [52] . 

In F-theory constructions one of the possible ways to break the GUT symmetries is 
to turn on a flux on the worldvolume of the seven-brane supporting the unified gauge 
group. In the present case, the SU{5) gauge symmetry can be broken by turning 
on a non-trivial flux along the hypercharge with Qy = diag{ — |, — ^j^j^}- As a 
result, SU{5) multiplets residing on certain curves where the fiux restricts non-trivially, 
might split. In the case of Higgs fiveplets in particular, this mechanism could be 
used to eliminate the unwanted triplets from the light spectrum, however the effective 
theory may no longer contain complete SU{5) representations and the proportionality 
conditions oc ka may not hold. We analyse this issue in the remaining of this section. 

Suppose that Miq, M5 are two integers representing the number of complete SU(5) 
10 and 5 representations in a specific construction. We work in the context of spectral 
cover approach, thus we consider the £"8 embedding of SU{5)gut and the breaking 

^8 ^ SU{5)gut X U{1)^ 

The SU{5) chiral and Higgs matter fields descend from the adjoint representation of 
the Es symmetry and reside on the various curves denoted with Sio^- , Sg, . We expect 
that the U{1) fluxes (those not included in SU{5)gut) together with the tracelessness 
condition ^ • — imply the following condition on the numbers of multiplets [T7] 

Y,M', + Y,Mio = (28) 

i j 

Consider flrst the case that we have all 10 type chiral matter accommodated only on 
one Eio curve and all chiral states 5 respectively on a single Sg curve. Then condition 
(f28]) imphes the relation Mio = -M5 = M. 

We denote with A^y^ , Ny^q the corresponding units of Y fiux which splits the SU (5) 
multiplets according to 



^_ . ( ™(3,l)_i/3 - "(3,1)1/3 = ^5 
^ ' \ "(1,2)1/2 - "(l,2)_i/2 =M5+ Ny, 



^^(3,2)1/6 - "(3,2)_i/6 - Mio 
"(3,l)-2/3 - "(3,1)2/3 = - Ny,, 

"(1,1)1 - "(i,i)-i = + ^y-io 

(29) 



Notice that these formulae count the number of 5-components minus those of 5 and 
the number of 10 components minus those of 10. Since we know that families are 
accommodated on 5's we require ^^^^^ > ?^(3,i)_j/3 which implies M5 < 0. Similarly, 
because the remaining pieces of fermion generations live on lO's, we wish to end up 
with 10-components after the symmetry breaking, hence we should have Miq > 0. For 
example, for exactly three generations we should demand Miq = — M5 = 3 and Ny^ = 0. 
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In general however, various curves belong to different homology classes and flux restricts 
non-trivially to some of them, thus A^y^. 7^ at least for some values of j. Furthermore, 
we note that in the above description we deal with negative numbers which are used 
to discriminate the 5, 10 conjugate representations from 5's and lO's. In order to take 
correctly into account their contributions to the /3-functions, care should be taken so 
that only positive numbers are introduced. Thus for the number of triplets we take 
I M5 1 = M, while for the required lepton doublets we take | M5 + Ny \ = M — Ny > 
and so onQ For this simple example we observe that the contributions to /3-functions 
come with the right proportionality factors 1 : 1 : ^ = : k2 : ki only in the unrealistic 
case Ny^ = Ny^g = N where |6y = 63 = 62 = 2M — ^ A^, so that the arbitrary cutoff 
A' can be removed to a redefinition of the gauge coupling as in the case of the gauge 
fields. 

In the general case the situation is more complicated. In the £"8 embedded SU (5) 
case in particular, in principle there exist ten S5 and five Sio curved where the cor- 
responding SU (5) representations may reside. The implementation of monodromy 
actions may reduce the number of these curves, however, the construction of realistic 
effective GUT theories capable of reproducing the know hierarchical spectrum eventu- 
ally require the involvement of matter fields arising from a variety of matter curves. 
Furthermore, because of ?7(l)y-flux effects several SU{5) representations split in a phe- 
nomenologically prescribed (and sometimes promiscuous) way resulting into incomplete 
multiplets. 

Bearing in mind the above general context, we write the contribution to the /3- 
functions as follows 



i j j 

62 = ^El^5 + A^5yl + ^El^i( 

i j 
&^ = ^El^5 + ^5>'l + ^El^! 



^ (30) 



+ I E l^iol + I E 1^10 - Ky\ + E 1^10 + ^^^oy I 
j 



where i,j summations are over the number of Sg and Sio discrete curves in the quotient 
theory, i.e. after the monodromy action. 

In order to acquire again the relation (1,1,5/3) we first must impose the equality 



^ In fact here we assume that IM5 1 counts exactly the number of (3, 1)1/3 s^nd that there are no extra 
(3, l)_i/3, (3, 1)1/3 pairs and similarly for the other components, otherwise threshold effects should be 
taken into account [30II31II33] . More general cases are considered in the next sections. 

^This is clear from the decomposition of the adjoint of Eg, under SU{5) x SU{5)± 

248 -> (24, 1) + (1, 24) + (10, 5) + (5, 10) + (5, TO) + (TO, 5) 
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between the contribution to SU (3) and SU (2) 

\ E l^5l + E l^iol + \ E = \ E l^'^5 + A^5y I + \ E l^/ol - ^ (31) 

where the parameter B was introduced for later convenience. We may rewrite the 
constraint above in a simpHfied way as fohows 

^ \Ml\ -Y.\Ml^ Nly\ = ^ \MU - ^ \M(, - Ni,y\ (32) 

i i j j 

Now, we can reexpress the U{1)y contribution in terms of B and a remaining quantity 
as follows 

by = ^B-2Y^ \M(,\ + \Mi, + Ni,y\ + ^ - AT^o^l (33) 
j j j 

In ()33p we have split the C/(l) /3-function into two parts, the first being ^B which 



comparing with (j3ip we observe that preserves the required fca-proportion with 62,3- 
Therefore, the remaining part must be zero and finally we get the following two con- 
straints 

^ \M(,\ - m = E i^i'o - ^iVi - E 1^5 + n5y\ 

^ ^ ^ ' (34) 

2 E i^ioi = E i^i'o + ^iVi + E 1^10 - ^iVi 

In writing the above contributions to the /3-functions we have assumed that they cor- 
respond to all the matter content, including the 3 generations of the SM. Clearly, it 
is a non-trivial task to obtain a low energy effective field theory model with only the 
MSSM spectrum. 

We note in passing that if we want to consider only the extra matter content, 
we should subtract a 3 from all absolute values: \Ml\ — > IM5I — 3, {M^ + N^y\ — )■ 
jMg + N^y\ — 3, etc. But since we know that the SM content forms a complete 5 and 
a 10, the contributions from the 3 generations obey the rule (1,1,5/3). Therefore, the 
above constraints hold also for the case when only the extra matter is considered. Note 
also that if all Ny- are zero the relations hold automatically, as they should. Of course, 
any choice of M's has to satisfy the constraint (p8]) . 



4 Application to a realistic minimal model 

From the analysis of the chiral and Higgs matter contributions to the /3-functions pre- 
sented above, it is clear that the cutoff independence constraints cannot be satisfied for 
any spectrum, even if flux and other consistency conditions are met. In this section we 
are going to discuss a realistic model with the minimal spectrum where these require- 
ments are fulfilled. We make an explicit calculation to determine the gauge unification 
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Table 1: Field representation content under SU{5) x [/(l)^. , their homology class and 
flux restrictions [T7] for the model [18] . Superscripts in the first column are numbering 
the curves, while subscripts indicate the family, the Higgs etc. For convenience, only 
the properties of 10, 5 are shown. 10, 5 are characterized by opposite values of U — )• —ti 
etc. Note that the fluxes satisfy N = Nj + Ng + Ng and Mio, + M5 . = while 

X = X7 + XS + X9- 



and b, t- Yukawa couplings through the renormalization group running in a model which 
fulfills all the requirements. 

To start with, we recall here the basic features of the model [18] which emerges 
from the general class [T7] presented in Table [TJ The first two columns give the field 
content under SU (5) x U (l)t^ for the case of Z2 monodromy. The third column presents 
the homology classes where ci is the first Chern class of the tangent bundle of GUT 
surface Squt and r/ = 6ci — t with —t being the first Chern class of the normal bundle 
to Sgut- The Xi are unspecified subject only to the condition x = X7 + X& + X9- If 
denotes the U{1)y flux, to avoid a Green-Schwarz mass for the corresponding gauge 
boson we must require Tyrj = J^yCi = 0. Then, we get Ni = J^yXi and consequently 

= J^Y ■ X = -^7 + -^8 + -^9- Using these facts, all remaining entries of column 4 in 
Tabled] are easily deduced. 

We now take the flux parameters to be Mio, 23 = !) Ms^ 24 = ~1 and = 0, while 
we have the freedom to choose Nf^s^g subject only to the constraint = + Ng + A'g. 
This choice of Mi,Nj^s ensures the existence of three 10 and three 5 representations 
which are needed to accommodate the three chiral families. 

Next we use the U{1)y flux mechanism to realise the doublet triplet splitting and 
make the model free from dangerous color triplets at scales below Mqut- We choose 
Ms^^ = 1, to accommodate the Higgs and Ms^^ = with A'g = 1 so that we get 
the splitting of down quark in the Higgs flveplet 

■/ ^(3.1)-i/3-^3,l),/3 =M55=0 

■ I ^(1,2),/, - ^(i,2)_,/2 =Ms,+iV-Ar8 = -l 
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In order to satisfy the trace conditions we choose Msg = — 1, A^7 = — 1 so that 5^^^ has 
only a colour triplet component: 



n 



-5(6) 



(3,1). 



1/3 



n 



(3,l)l/3 



"(l,2)l/2 -"(l,2)_i/2 



= -1 

Msg + iV - iV7 = 



(36) 



In this simple example we have succeeded to disentangle the triplet from the up-Higgs 
curve at the price of generating however a new one in a different matter curve. Yet, 
this allows the possibility of realising the doublet-triplet splitting since we can gen- 
erate a mass Mt for the triplet by coupling it to an antitriplet via the appropriate 
superpotential term [H]. This way we obtain the corresponding Higgs doublets light. 

However from Table [J one may see that the matter on the S^q{2,3) curves will be 
affected by the A^y^g flux. In particular the content of 10/10-representations on S]^q(2,3) 
is split by the choice of fluxes as follows 



^10(2) : ' 



n 



(3,2)i/6 - "(3,2)_i/6 

"(3,l)_2/3 - "(3.1)2/3 



I "(1,1)1 



n 



(1,1)- 



: < 



n 



(3,2)i/6 



n 



n 



(3,1)- 

(l,l)l 



2/3 



(3,2)_i/6 

"(3,1)2 



^2/3 



n 



(1,1)- 



: M1O2 = 1 

= M102 + ^7 

Mio, = 1 

M103 - ^8 ■ 

Mio, + Ns - 



■- 2 

: 


2. 



(37) 



(38) 



We observe that in the presence of flux one = (1, l)i state is 'displaced' from S 



to the S]^Q(3) curve. A similar dislocation occurs for one 



(3,1) 



-2/3 



of S 



10(3) 



10(2) 

which 



'reappears' in ^iq(2) ■ We note that this fact implies a different structure for the up, down 
and charged lepton mass matrices. It can be checked that the particular distribution 
of the chiral matter on the specific matter curves can lead to interesting results with 
respect to the fermion mass structure and other phenomenological properties of the 
model [18]. For clarity, the final distribution of the MSSM spectrum along the available 
matter curves is summarized in Table EJ 

Here, we are interested for the renormalization group properties of the model. Using 
the above context we get the following relations 



i i 

j j 3 

which obviously respect the three constraints in (I34p and at the same time 



(39) 



(40) 



which ensures three chiral families and the tracelessness condition (|28 



13 





M 


N 


Q 








M 


N 


dP/T 


r 1 'U 


10(1) (F3) 


1 





1 


1 


1 


5(0) 


{hu,T) 


1 





1 


1 


10(2) (F2,i) 


1 


-1 


1 


2 





5(1) 


ih) 


-1 





-1 


-1 


10(3) (Fi,2) 


1 


1 


1 





2 


5(2) 


ih) 


-1 





-1 


-1 


lo(^) (-) 

















5(3) 


(-) 


























5(4) 


(/i) 


-1 





-1 


-1 














5(5) 


(hd) 





-1 





-1 














5(6) 


(T) 


-1 


1 


-1 






Table 2: The distribution of the chiral and Higgs matter content of the minimal model 
along the available curves, after the U{1)y flux is turned on. The three families Fi = 
10i,/j = 5j are assigned on the curves as indicated. The Higgs doublets hu,d and T/T 
triplets descend from three different curves, (see also Table [J and text.) 

5 The GUT scale and the flux thresholds 

In the previous section we have shown that the contribution to the /3-functions of the 
chiral spectrum of the particular SU{5) model is in the required ratio, i.e. proportional 
to the coefficients ka- 

We have already explained that in the present model we assume that the breaking 
of SU{5) occurs when a flux is turned on along the U{1)y component of the SU{5) 
gauge symmetry. The rather interesting feature of the particular choice of fluxes leads 
to a minimal spectrum which consists of the MSSM content and only a color triplet 
pair below Mqut- Furthermore, the appearance solely of (any number of) pairs of 
extra triplets -and no other additional matter- has the interesting property of not 
disturbing the scale where the uniflcation is achieved. This fact is valid irrespectively 
of the scale that these triplets become massive. As a result, the nice properties of 
the minimal supersymmetric unification scenarios where the GUT scale is determined 
around Mqut ~ lO^^GeV is retained. On the other hand, the value of the gauge 
coupling gcuT does depend on the number of triplets and the scale they become massive. 

For convenience we start by recapitulating the basic analysis of the gauge coupling 
running in the presence of fluxes. In F-theory constructions the flux mechanism em- 
ployed to break the GUT symmetry induces a splitting of the gauge couplings at the 
GUT scale [5l[29]. The following relation is found at the unification scale |29j 

1 _ 5 1 _ 1 ^2 1 ^^^^ 

aviMGUT) 3 aiiMcuT) a2{MGUT) SasiMcuT) 

The running of the couplings is governed by the RGEs where contributions to the /3- 
functions come from the MSSM spectrum and other possible remnants from the higher 
theory. In the case of the model under consideration we need to consider the simple 
case where in addition to low energy states only an additional triplet pair appears in 
the spectrum below the 5[/(5)-GUT breaking scale Mqut- We assume that at some 
scale Mx < Mqut the extra triplet pair T, T decouples and only the MSSM spectrum 
remains massless for scales fx < Mx ■ The low energy values of the gauge couplings are 
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then given by the evolution equations 



' +^ln^ + ^ln^ (42) 



aa{Mz) aaiMcuT) 2tt Mx 2tt Mz 

where 6^', {ba) are the /3-functions above (below) the scale Mx- 

Using the GUT relation (j4ip one arrives at the following formula for the GUT scale 

2i„ fMxV'^ 



Mgut = ef^^" i^-j^ j Mz (43) 

where ^ is a function of the experimentally known low energy values of the SM gauge 
coupling constants 

' ' ' 12 1 ^^^^ 



A 3ai(Mz) a2{Mz) SasiMz] 



and p is the ratio 



"^1 <^^> 

where /3, are the ^-functions combinations above and below the Mx scale 

/3. = 6f--6f-^6^ (46) 

/3 = bY-b2- (47) 

If the only additional states above Mx are the color triplets T = (3, 1)1/3, T = 
(3, l)_i/3, it can be easily checked that their total contribution to (3x combination 
adds up to zero. This means that fix = fi and p = 1, thus Mqut in (|^3]) becomes 
independent of the Mx scale and in fact it is identified with the MSSM unification 
scale 

Mu = Mgut = Mz ~ 2.15 x 10^'^GeV 

This result is only slightly modified when two loop corrections are taken into ac- 
count [30]. The irrelevance of the triplets' decoupling scale Mx holds for the one 
loop calculation and it is adequate for our purposes. With these preliminaries we are 
ready to discuss the computation of the Yukawa couplings. 



6 Wavefunction overlapping integrals 

In the class of models under consideration U{1) symmetries are acting as family symme- 
tries and fermion families are distributed over the curves in such a way so that fermion 
mass textures are rank one. Then the only tree level Yukawa couplings are those of 
the third generation and these are computed in terms of integrals over the overlapping 
wavefunctions at the triple intersection points 

Xij = M^ / ipi'ifj(f) dzi A dzi A dz2 A dz2 (48) 
Js 
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Chiral matter is localized along the matter curves and therefore the value of the integral 
involving the trilinear (top and bottom) Yukawa couplings principally depends on the 
local shape of the corresponding wavefunctions. In addition, a crucial factor in the 
determination of the exact values is the global normalization of the wavefunctions 
which could in principle depend on the geometry of the compact manifold. In what 
follows, we will make use of the analysis of the previous section to argue that this 
dependence enters only through the torsion which according to Ray-Singer theorem is 
independent of the metric of the manifold. 

To compute the integral we use the knowledge of the wavefunctions close to the 
intersection point. The solution of the zero-mode equations leads to a Gaussian profile 
of the wavefunctions ip which for localized solutions acquire the general form 

oc e""" ^ (49) 



where m is a mass parameter related to some background Higgs vev and q = \/ q'f + g^- 
We also need the normalization of the wavefunction which is determined by the integral: 

c = Mt / l^pdz ^dz = TT^n'^ (50) 

Js ^ Q 

The factor vr— 3^ is the result of the gaussian integration along the coordinate normal 
to the curve. The factor Ti? is introduced to account for the integration along the 
coordinate parametrising the curve. We observe that in principle three different scales 
are introduced in the above normalization formula of the wavefunction, namely M*,m 
and TZ^^. Thus, although the solution of the equations of motion lead to a wavefunc- 
tion which is peaked on the curve while falling off exponentially away from it, the final 
profile appears to comprise global information through the aforementioned mass pa- 
rameters. We argue here that these apparently unrelated mass scales are connected by 
our previous analysis and the final wavefunction formula is solely characterised from an 
invariant quantity related to the torsion. Indeed, the scale m that has been introduced 
in relation to the vacuum expectation value of a background scalar field in connection 
with the breaking of the enhanced gauge symmetry, is associated to the scale M^, ~ Mq 
which appeared in ()24p . On the other hand, the parameter TZ 'measures' the integration 
along the matter curve inside the GUT surface S, thus clearly Tl~^ < Mc and from our 
previous analysis one naturally expects that Tl~^ ~ Mkk ~ Mgut- Although both 
scales are not known exactly, it is possible from our previous analysis on threshold 
corrections to obtain a scale independent normalization of the wavefunction. Indeed, 
following the reasoning presented in the above lines we can write (j50p as follows 

C = ^4^ (51) 



Prom the Mgut definition formula (j23p we observe that the ratio of the two scales 
Mgut/Mc is independent of the geometry of the manifold and can be written as a 
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simple exponential function of the torsion. Substituting ()24p into (j5ip we find that the 
wavefunctions must be normalised with 

J- = , /i e^/SC^s/e-To) 
VC Vvr 

Therefore, the normalization constant C is found to be independent of the two scales 
Mgut, Mc of the theory. 

Now, the computation of the integral for three generic wavefunctions of the form 
()48p with arbitrary charges qi,q[,q'l,i = 1,2 participating in the triple intersection the 
Yukawa couplings, gives 

q + q' + q" {qiq'2-q'iq2? ^ ' 

where q'f = —qi — q[ from charge conservation in the triple intersection, while q = 
\/ qf + q'2 and similarly for q',q" [18]. We point out that the dependence on the value 
of gauge coupling acuT in previous wavefunction normalizations has been replaced 
by the torsion. Previous normalizations where plagued by the smallness of the gauge 
coupling constant acuT ~ ^ which led to rather small values of the top Yukawa 
coupling in disagreement with the data 0. The hope is that the replacement of this 
small factor with the exponential factor incorporating the analytic torsion might point 
to the correct answer. 

In order to demonstrate the points discussed above, we proceed with the compu- 
tation of the determinants specifying the value of the torsion in a simple case. As in 
ref [5] we take a line bundle 0{n, —n) on a Hirzebruch surface Fq = x P^. Taking 
into account that the Euler character on P^ is given x(-P^) C'(ra)) = n + 1 the product 
formula of Ray-Singer torsion yields 

To{n,m) = in+ l)To[m) + {m + 1)7^(„) 

Since we want to eliminate the color triplets (3,2)_5/g + (8,2)5/5 we need to take 
xiP^, C'(n)) = n + 1 = 0, therefore for the line bundle C^^^ we take n = — 1, so 

"Ts/e = "70(1,-1) = 27c:)(_i) + (-1 + l)7ci(i) = 2To{-i) 
For the trivial line bundle we take n = 0, hence 

To = 70(0,0) = 27c)(o) 

We now recall that a positive elliptic differential operator, as is the case of 7^ 
with spectrum A„, is associated to the zeta functions [MjlSS] 



Ckis) = 



7n 

As 



^ In previous normalizations, the appearance of the gauge couphng in the wavefunction normalization 
is due to the relation a^J^j, — Mf J cj A cj oc TZ'^Mf. However, in this approach, the result depends on 
the scaling of the metric and it is not clear whether the effective value gauT or the coupling gs should 
be involved. 
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where 7„ is the degeneracy of A„. The determinant is defined by 

dCk{s) 



log(DetA) 



ds 

1/-' I 



s=0 



Therefore, the torsion is given 7o(fc) = — ^Cfe(O) with 

fc+i 



cm = <r{-1) -l{k + l) + ^(2/ -\k + l\) log / 

1=1 

where Ck(-1) ~ -0.165421. Then we get 

Ts/e - To = 2(ro(_i) - Toio)) = Co(0) - CM = -\ 
For this particular choice of bundle, the Yukawa coupling then becomes 



, 40F/e {qq'q"Yl^ 
q + q' + q" {qiq'^ - q[q2Y 



6.1 The top mass 



Let us now apply the above formula in the case of the model presented in the previous 
section. From Tabled] we see that a tree-level mass for the top quark is available from 
the gauge invariant Yukawa coupling 

t^j^ee = Atl03 • IO3 • hh^ 

ti t2 — tl — t2 

where from the second line we see that the U{1) invariance Yli^i = is also satisfied. 
The charges in (j53p for the top vertex for this model have been calculated in |18J 




M''?^} = J^, -j= , {q[,q'2} = \ T^,-^ (54) 




The resulting top quark coupling is computed form ([53]) and is found to be At ^ 
1.23. Similarly the bottom quark Yukawa coupling in the model is obtained form the 
superpotential term 

Wtree = Afe IO3 • 83 • 

t2 h + ^5 ^1 + *4 

while repeating the above analysis one finds A^ ~ 1.17. We thus observe that we are in 
the large tan /3 regime with a top coupling in the range Aj > 1 which implies a physical 
mass around 186 GeV which is within the experimental range |36j . 

In deriving the above numerical values we have made a specific simple choice of 
flat line bundle to compute the torsion which is involved in the top and bottom mass 
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0.6 0.8 1.0 1.2 1.4 1.6 1.8 



Figure 1: Plot of the top mass versus the value of the Yukawa coupling Xt{MGUT), for 
tan/3 = (40 — 55). Horizontal lines indicate the experimental uncertainties of the top 
mass. 

formula through the exponential factor appearing in ()52p . Of course, one may think of 
more involved cases of bundle structures however the above estimates are not expected 
to alter dramatically. Moreover, in a fully realistic model other parameters emerging 
from various thresholds in the renormalization group running etc may modify also the 
computation in obtaining the low energy masses. The final result will also depend on 
the specific choice of tan /3 which is fixed only by phenomenology and not by some the- 
oretical principle. Thus, in our present computations we have the freedom to adjust its 
value to fit the low energy mass data. Furthermore, the reliability of our computations 
is corroborated from the rather interesting property of the top coupling, namely its 
fixed point behavior. This means that for a range of relatively large initial values 
at GUT, the renormalized value of the top Yukawa approaches the same limiting value 
at the scale of M^y. It is anticipated that another choice of bundle would imply rea- 




Figure 2: Plot of the bottom mass versus the value of the Yukawa coupling Xij{Mgut) 
for three different values of tan/3 = (40,45,50). Horizontal lines indicate the experi- 
mental uncertainties on the bottom mass. 
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sonable deviations from the numbers obtained in this particular example, so that they 
do not alter substantially the low energy predictions. In this sense, a fully realistic 
model could arise under an appropriate type of bundle structure. For completeness, 
we demonstrate the fixed point behavior by a one-loop renormalization group analysis 
using analytic formulae [39] of the top-bottom Yukawa coupled differential equations 
in the Appendix. To visualize this effect, using the derived formulae, in Fig.([T]) and 
([2]) we plot the low energy top and bottom masses as a function of the values of the 
corresponding Yukawa couplings at Mqut- This is done for a range of large values of 
tan /3 w (40 — 55) as indicated in the plots. We can see clearly in the bottom mass 
plot the three curves corresponding to the values tan j3 = 40, 45, 50. Notice that in the 
case of the top mass plot, all tan/3 curves are squeezed to one because the top mass 
depends on sin(/3) < 0{1). Horizontal lines specify the experimental uncertainties. It 
can be seen that top mass is found in a region where the Yukawa coupling is fairly 
close enough to the fixed point and small uncertainties in its computation are within 
the experimentally acceptable bounds. 

7 Conclusions 

In this work we have reconsidered the possibility of determining in a reliable way the 
mass scales and moduli parameters which enter in the computation of the tree level 
Yukawa couplings in a class of F-Theory GUT models. We consider in particular 
the SU{5) GUT model generated by a seven brane wrapping a complex surface S 
and embedded in the exceptional Es underlying theory. The SU{5) GUT breaking 
is supposed to occur when a U{1)y flux is turned on. In these models chiral fields 
are accommodated along the Riemann surfaces which constitute the intersections of 
the GUT brane with other 7 branes where the SU{5) gauge symmetry is enhanced. 
Chiral matter fields lying on these 'curves' are charged under additional C/(l)'s which 
act as family symmetries. Under certain monodromy actions, Yukawa couplings can 
be formed at triple intersections of such curves and predict rank one fermion mass 
textures. The values of the Yukawa couplings are expressed in terms of integrals of 
overlapping wavefunctions of the matter fields participating in the triple intersection. 
The wavefunctions are found from the solution of the equations of motion and their 
normalization in principle appears to require knowledge of the mass scales of the theory 
and topological properties of the compact space and matter curves. We have made 
reasonable assumptions about the mass scales of the theory and used previous analysis 
of threshold effects from KK modes [Slj^^T] to determine the GUT scale and the mass 
ratios which control the wavefunction normalization. Further, we have studied similar 
effects from a generic chiral sector of the F-theory model in the presence of C/(l) fiuxes 
and have analysed the constraints imposed on the light spectrum of the effective field 
theory model. We have incorporated the above findings into the computation of generic 
trilinear top and bottom Yukawa couplings in the context of the SU (5) model embedded 
in Es and found that the couplings depend only on the extra U{1) charges carried by the 
fields at the intersection and a topological invariant quantity, namely the Ray-Singer 
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torsion [28]. Furthermore, we have apphed the method of our computation in a viable 
model obtained in the described F-theory context [18J with the minimal light spectrum 
which meets all the requirements imposed by our analysis. To compute the torsion, we 
have used [5] as example a simple case of flat line bundle on a Hirzebruch surface. We 
have performed the renormalization group analysis, taking into account the KK-mode 
effects discussed previously as well as flux thresholds to determine the low energy values 
of the gauge and Yukawa couplings. We have found a top and bottom mass which is 
in agreement with their experimentally determined values. 
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Appendix 

In subsection 6.1 we applied our method for the computation of the top and bottom 
Yukawa couplings which for a specific choice of line bundle they where found to obtain 
0(1) and comparable values (large tan/3 regime) at the GUT scale. The top and 
bottom masses are given from their renormalised values at scales ~ Miy, and are 
obtained from a solution of two coupled differential equations (DEs) (we assume for 
simplicity small tau Yukawa coupling). For comparable initial Xt^b values the DEs 
cannot be disentangled, thus in order to give a reliable calculation of their low energy 
values, we need to solve the above DEs and plot the values of the top and bottom 
masses for reasonable ranges of X^^t initial values. To this end, in this Appendix we 
obtain analytic formulae by extending the analytic solution [37] of the system of the 
1-loop RGE for the top and bottom Yukawa couplings ht and h^. Defining at = /(47r) 
and ab = A^/(47r), the RGEs are 



= ^ (6at + - Gq) at 

^ = 77- ("t + - Gd) Q-b 
at 2tt 



(55) 



with Gj = c{o;i + C2a2 + c^a^, I = Q,D, and the constants c are given by {SU{5) 
normalization) 



,1^2 , 



} 



16 13 
3'"^' 15 



and {cf,c^,cf} 



1^ q 1 
3 ' ' 15 



Ignoring the small difference in the U{1) factor of the constants c, we shall make the 
transformations 



at = JqX, 



where 



7q = exp 



and ab = ~ -fqy 



^ fGQ{t')dt' 
27r J to 



(56) 



The system (j55|) becomes 



Simple manipulations give 



dx 


1 


IE 


~ 2^ 


dy 


1 


dt 


~ 2^ 



(6x + y) 7qX 
{x + 6y) 



(57) 



Therefore 



7 



dln{x — y) 
~dt 



^^^%^ = ^6(x + y) 
dt 2-K ^ ^' 

dln(xy) 7q 

d\Yi{xy) f xy \^ 

dt ' \xoyoJ 



6 



x-y 
xo - yo 



(58) 



(59) 
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where xq and uq are initial conditions at the scale to- Observe that, since 7Q(to) = 1) 
we have that xq = at (to) and yo = ab{to). Using ([59]) we can write 



(x + yf = {x- yf + Axy = [x - yf + 4. 



xoVo 



and the first of the equation in (l58l) reads 



d{x -y) _1q 



xoyo 



Defining a further new transformation 

u = k^l^^l^ = ko/{x - yfl\ 



where h 



■o 



Xoyo 



'{xo-yoy/^ 



we get 



du 
It 



27rwii/6 



27r 



12/5 



+ ko 



u J 



7/5 



5 2 7,6/5 VTTu 



,2 r,6/5 



ni/5 



Therefore 



,1/5 



Integrating the left-hand side we get 

U^^^ ^ 10 1/5 

du= — u^'^ 
7 



du = -^kl^^Jqdt 



1 1 6 



(60) 



(61) 



where 2-^1 is the hyper geometric function. Let us work on the right-hand side of (I6ip 



iQ = exp 
Since 

we get 



ai{t')dt' 



to 



exp 



^ 1^ (c?ai(t') + c^a2(t') + c^a3(t')) 



to l-^(t'-to) 



27r 

-—In 



27r 



-(t-to) 



1=1 



l-^(t-to) 
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Define for convenience tlie new constants (we suppress liereafter the superscript Q) 

biaio/{2ir) = Bi, and Ci/bi = ei 

Integrating 7q we get 
1 



rgdt 
Bi 



Bi(l + ei] 



-1 + Bi{t - to)] [1 - Bi{t - to)r [1 - B2{t - to)r X 



{l-B2{t-to)) 



Bi — B2 

[1 + ei, -62, -63,2 + ei 



[l-B^{t-tQ)] 



63 



Bi 



{l-B3{t-to)) 



Bi — Bs 

B2{-1 + Bi{t - to)) B^i-l + B^{t - to)) 



-63 



Bi — B2 



Bi — B3 



where aJ' is the Appell hypergeometric function of two variables [38]. Simphfying the 
relation and using 1 — Bi{t — to) = aio/ai{t), we get 



ai{t) 



ei+l 



aJ^ { l + ei,-e2,-e3,2 + ei,' 



Bi — B2 


62 


Bi - B3 


[ Bi \ 




[ Bi \ 


"10 ^ 


Bz 


( "10 



63 



B2 - Bi \ai{t) J' Bs-Bi \ai{t) 
Therefore the differential equation (j6ip gives 
10 



7 



■ u 



1/5 



, ^ / 1 1 6 

Vl + u- 2-7^1 ( 5' 2' 5'~^ 



1 



2tt ° Si(l + ei) 



"10 



ai{t) 



61+1 



Bi — B2 


62 


Bi — B3 


63 


[ Bi \ 




[ Bi \ 





aJ^ { 1 + ei, -62, -63, 2 + ei. 



Bo 



"10 



S3 



52-^1 \ai{t)J ' B^-Bi \ai{t 



aio 



+ constant 



and simplifying 



/ ^,116 

Vl + n - 2-Ai 7, -, 7, -u 



5' 2' 5' 



7 ,6/5 1 



47r ° Si(l + ei] 



QlO 
ai(t) 



61+1 





Bi — B2 


62 


Bi — B3 




[ Bi \ 




[ Bi \ 



63 



(62) 



J" 1 + ei, -62, -63,2 + ei. 



B2 ( aio 



S2 - Bi Vai(i) J ' B3-B1 \ai{t) 



+ constant 



Let us see now how x and y are related with a; which, from (j60p is given in terms of ?/. 
Since (see (fSQil ) 

/ a; \ '^/^ 

X - y = X + (-y) = (J, and (-y) • x = -xoyo ( — j 
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we get (using ([60|) for the second equality) 



2 



x = - [1 + ^/l + kQuF^ 



-1 + Vl + A;ocj-5/6 



2\u 
2\u 



j (1 + 



6/5 



-i + ^/TT^) 



(63) 



From ([56|) we get 



= 7q a: 



ab = lDy 



n 

3 

n 



l-^(t-io) 



27r 
2tt 



{t - to) 



1 /A;o 



2 V u 



6/5 



(64) 



These formulae incorporate the RG running effects and encompass all the information 
for the initial values of the gauge and Yukawa couplings at the GUT scale. We may 
substitute the top and bottom initial values and use them ^9] to derive our results. 
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